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ABSTRACT 

Virial mass is used as an estimator for the mass of a dark matter halo. However, the 
commonly used constant overdensity criterion does not reflect the dynamical structure 
of haloes. Here we analyze dark matter cosmological simulations in order to obtain 
properties of haloes of different masses focusing on the size of the region with zero 
mean radial velocity. Dark matter inside this region is stationary, and thus the mass 
of this region is a much better approximation for the virial mass. We call this mass 
the static mass to distinguish from the commonly used constant overdensity mass. We 
also study the relation of this static mass with the traditional virial mass, and we find 
that the matter inside galaxy-size haloes (Af « IO^^Mq) is underestimated by the 
virial mass by nearly a factor of two. At z « the virial mass is close to the static 
mass for cluster-size haloes (M « IO^^Mq). The same pattern - large haloes having 
Mvir > Mstatic " cxists at all redshifts, but the transition mass Mq = Mvir — Mstatic 
decreases dramatically with increasing redshift: Mo{z) « 3 x 10^^/i~^Mq(1 -t- z)~^-^. 
When rescaled to the same Mq haloes clearly demonstrate a self-similar behaviour, 
which in a statistical sense gives a relation between the static and virial mass. To our 
surprise we find that the abundance of haloes with a given static mass, i.e. the static 
mass function, is very accurately fitted by the Press & Schechter approximation at 
z = 0, but this approximation breaks at higher redshifts z ~ 1. Instead, the virial mass 
function is well fitted as usual by the Sheth & Tormen approximation even at z < 2. 
We find an explanation why the static radius can be 2-3 times larger as compared with 
the constant overdensity estimate. The traditional estimate is based on the top-hat 
model, which assumes a constant density and no rms velocities for the matter before 
it collapses into a halo. Those assumptions fail for small haloes, which find themselves 
in environment, where density is falling off well outside the virial radius and random 
velocities grow due to other haloes. Applying the non-stationary Jeans equation we 
find that the role of the pressure gradients is significantly larger for small haloes. At 
some moment it gets too large and stops the accretion. 

Key words: dark matter - galaxies: haloes - large-scale structure of Universe - 
cosmology: theory - methods: iV-body simulations 



1 INTRODUCTION 

The c urrently accepted p a radigm of hierarchical clus- 
tering (jWhite fc ReesI 1 19781 . iBlumenthal at al.l 1 1984 ) pro- 
vides a picture of assembly of dark matter haloes in 
which more massive haloes are formed through merg- 
ing and accretion of smaller ones. This picture is sup- 
ported for example b y recent observations of o n going merg- 
ers i n clusters (e.g. Ivan Dokkum et al. 1 1 19991 . iRines et al.l 
l2007f) . The theoretical framework of the cold dark mat- 
ter (CDM) has been proved successful in the descrip- 
tion of the structure formation in the Universe from tiny 
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fluctu ations in the primordial density field (see IPrimaclj 
l2003l for a review). This model has received support from 
many diff'erent obs ervations, in part i cular of the gravita- 



tional lensing effect llSmith et al ]|200ll.lGuzik fc Seliay 
iKneib et al.]|2003l,lHoekstra et al.ll2004 ISheldon" et al.ll 



2002 . 



2004 . 



iMandelba um et al.|[2006l). CMB ( ISpergel et al."200 7|), the 



abundance of clusters jPierpaoli et al. 2001, Gladdc rs et al.l 
20071). and satel lite dynamics fZaritskv fc Whit JTToM 



Prada et al.ll200^ ). Cosmological simulations with ever in- 



creasing resolution play important role by making accurate 
predictions of different properties of dark matter haloes. 
Results of those simulations are used by other methods. 
For example, semi-analytical models of galaxy formation 
have been either incorporated into N-body simulations or 
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use statistics such as halo mass function and merging trees, 
which were calibrated and tested using the simu lations (e.g. 
ISomerville fc Primackl Il999l : ICroton et al.l l2006h . The sim- 
ulations reveal important information about the internal 
structure of dark matter haloes (e.g . Navarro et al.l Il997l . 



iBullock et al]l200ll . IXavlor fc Navarrdl200ll ') which is reflect- 
ing their underlying dynamics. 

Many results of simulations implicitly use some defini- 
tion of what is the size and mass of a collapsed dark matter 
halo. The problem is that there is no well-defined bound- 
ary of a halo: density field is smooth around the halo. The 
common prescription for this boundary (and hence the mass 
belonging to the halo) is defined through the spherical col- 
lapse model (|Gunn fc Gottlll972l . lGunnlll977l ). The size of 
the halo at redshift z is given by the half of the radius of 
the spherical shell at turnaround which is collapsing at that 
redshift. This is the virial radius. Thus, it is very common 
to measure the mass of haloes in cosmological simulations 
taking the particles inside a sphere of fixed spherical over- 
density or to take those particles which are connected by 
a inter-particle separation below a given value (the friends- 
of-friends algorithm). We note that there is little justifica- 
tion for using the top-hat collapse model. Haloes do not col- 
lapse from perfect spherical homogeneous distribution. The 
environment of haloes is typically very non-spherical with 
most of accretion happening from few elongated filaments. 
The random velocities of the accreted matter also cannot 
be neglected. As the fluctuations collapse, the dark matter, 
which is being accreted, increases it rms velocities. This ef- 
fective pressure should affect the accretion rate. The only 
motivation for using the top-hat model comes from simula- 
tions. Indeed, early simulations indicated tha t the radius of 
overd ensity 178 is close to the virial radius (ICole fc LacevI 
1 19961 ): " the radius rirs approximately demarcates the in- 
ner regions of haloes at r < rirs which are in approximate 
dynamical equilibrium from the outer regions at r > ri78 
which are still infalling". Thus, the radius of overdensity 
200 (« 178) became the virial radius. The reason why this 
was a good approximation is simply coincidental: the early 
simulations were mostly done for cluster-size haloes and, in- 
deed, for those masses the virial radius is close to the radius 
of overdensity 200. The early models were fiat models with- 
out the cosmological constant. Models with the cosmological 
constant have produced significant confusion in the commu- 
nity. The top-hat model must be modified to incorporate the 
changes due to the different rate of expansio n and due to th e 
different rate of growth of perturbations fsee |Primacklll997[ ). 
That path produced the so called virial radius, which for the 
standard cosmological model giv es the radius of overd ensity 
relative to matter of about 340 (|Bullock et al.ll2001^ . StiU, 
a large group of cosmologists uses the old overdensity 200 
relative to the critical density even for the models with the 
cosmological constant. 

In this paper, we cast some light on this subject by 
searching for a physical extent of dark matter haloes, which 
is related to the physical processes t hat occur around col- 
lapsed structures. iPrada et al.l (|2006l ) provided first results, 
which indicated that spherically averaged, mean radial ve- 
locity profiles show an inner region in which there is no net 
infall or outfiow. The size of this region in virial units is 
mass-dependent: for galactic-size haloes it may even reach 
three times the virial radius. We use this result to study 



the properties of the mass inside this region and in particu- 
lar, to determine the redshift evolution of such mass. Recent 
effort (|Wechsler et al.ll2002l . fzhao et al.ll2003l ) has been de- 
voted to the analysis of the evolution of virial mass inside 
dark matter haloes, sometimes referred to as the mass accre- 
tion history. However, the picture presented in these works 
has been recently put into question. The analysis of the 
mass inside a fixed physical radius has revealed that galaxy- 
size haloes experience unphysical growth of their virial mass 
jPiemand et al.ll2007^ . The mean background density de- 
creases as the Universe expands, making the virial radius to 
increase even in the case of no accretion or mergers. This 
is clearly an artifact of the definition. Other works realized 
this issue, for example, the evolution of the spin parame- 
ter when the matter inside a fixed radius is taken into ac- 
count differs from that when using the evolving _Rvir instead 
(D'Onghia & Navarro 2007). On the other hand, there is an- 
other effect apart from accretion and merging, which has not 
been taken into account: haloes may also grow via relaxation 
of the surrounding regions near them. This is an important 
effect as it is a refiection of the dynamical processes which 
are turning non-virialized mass in the outskirts of a halo 
into the mass associated to it. Moreover, the long-term evo- 
lution of dark matter haloes show an interesting feature for 
ACDM cosmologies: their mass turns out to converge to an 
asymptotic val ue which depends on the definition of mass, as 
pointed out bv lBusha et al. (2005). Thus, it is desirable that 
the mass of a halo is measure d using a virialization-based 
criteria (see lMaccio et al"]|2003l for an interesting approach), 
instead of using boundaries of a given overdensity. 

The interest of the measurement of the physical mass 
associated to dark matter haloes is not only theoretical. In- 
deed, it may have a great impact on the number of col- 
lapse d objects in a given range of mass, i.e. the mass func- 
tion (|Whitell200ll ). Besides, it is of great importance for the 
process of formation and evolution of galaxies, and, hence, 
it is relevant for the results obtain ed from semi-analy tical 
modelling of galaxy formation (e.g. ICroton et al]|2006l ). to 
predict the main properties of observed galaxies. Thus, it 
turns out to be mandatory to bring attention on the mass 
belonging to a halo, if accurate predictions of the physics be- 
hind galaxies from cosmological simulations are to be drawn. 

This paper is organized as follows: in Section [5] we de- 
scribe the set of cosmological simulations used in our anal- 
ysis and the properties of the halo samples. In Section [3] we 
present the definition of static mass and how it is related 
to objects of different sizes. A brief analysis of equilibrium 
in this context is presented in Section m The scaling prop- 
erties of the static to virial mass relation with redshift are 
shown in Section (5] together with a simple model derived 
from this scaling relation. We obtain the static mass func- 
tion and compare it with analytical models in Section [S] In 
Section [7| we present the evolution of the static mass track- 
ing the halo progenitors. We discuss our results in Section [S] 
and present our conclusions in Section |51 



2 iV-BODY SIMULATIONS AND HALO 
SELECTION 

We use four different high-resolution ACDM simulations, 
which were run using the Adaptive Refinement Tree (ART) 
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code l|Kravtsov et al.lll997l ). These cosmological simulations 
are selected in order to span a broad range of halo masses 
with very good statistics of the number of haloes. The values 
of the cosmological parameters are flo — 0.3, fJbar ~ 0.045, 
JIa ~ 0.7, as — 0.9 and h — 0.7. Main simulation parameters 
may be found in Table [T] 

The resolution and the halo mass range covered by these 
simulations allow us to perform a detailed analysis of the 
properties of dark matter haloes with very good statistics. 
Dark matter haloes are ide ntified using the Bo und Density 
Maxima halo finder fBDM. lKTvpin et al.lll999l '). In order to 
resolve the inner structure and dynamics of dark matter 
haloes, only haloes which enclose more than 2, 000 particles 
inside the virial radius at 2 = 0, were chosen. We take into 
account all the mass (both bound and unbound particles) 
inside the virial radius of the halo. This radius is defined 
here as the radius of the s phere enclosing an overdensity 
given bv lBrvan fc NormanI ()l998h . i.e.: 

Ac = IStt^ + 82a: - 39x^ x = n{z) ~ 1 (1) 

In Figure [T] we plot the virial mass function of dark 
matter haloes over five orders of magnitude. Here the er- 
ror bars represent the Poissonian error, showing very good 
statistics for all except for the highest mass bins. The differ- 
ent symbols shows which simulation is used for each mass 
bin. The relation between mass bins and simulations de- 
pends on the fact that we only choose haloes enclosing more 
than 2, 000 particles inside the virial radius at 2 = 0. If two 
different simulations pass this constraint, then we choose 
the one with hig hest statistics. As previously shown (e.g. 
ISheth et al. I I2OOII ). the virial mass function in cosmological 
simulation s is well fitted by the Sh eth and Tormen function if 
a ~ 0.707 (|Sheth fc Tormenll 19991 . hereafter ST) in the mass 
range under study. These four simulations, as mentioned 
above, cover a broad range in mass, from lO^°/i"^M0 low 
mass haloes to cluster-size haloes with nearly 10^^/i~^Mq. 

In order to select a sample of haloes for our analysis, 
we take all the haloes whose centres are not inside the virial 
radius of a larger halo (hereafter distinct haloes). This crite- 
rion reduces strong environmental effects: internal dynamics 
of haloes is less affected by massive neighbours. 



3 STATIC MASS OF DARK MATTER HALOES 

In this Section we introduce and motivate our definition of 
static mass. When using cosmological simulations, it is com- 
mon to define the mass belonging to dark matter haloes in 
terms of a given spherical overdensity or inter-particle dis- 
tance. However, so far, it is computationally challenging to 
define a real virial mass, i.e. the mass which is showing fea- 
tures of statistical equilibrium, for a large number of haloes. 
This has to be kept in mind, if accurate results are to be 
found, although there is some kind of agreement between 
different definitions based on overdensity (| White! [20011 ). A 
particular case in which an overdensity-based definition may 
fail in providing a faithful description of dark matter haloes 
is the mass accretion history. Figure [2] clearly demonstrates 
this. Here we show the evolution of the density profile of 
a very high resolution (more than two million particles in- 
side the virial radius) galaxy-size halo from our cosmological 
simulation Box20, which has a better spatial (0.152 h~^kpc) 




10-' I 

10'° 10" 10'= 10'= 10'" 10'5 

Virial Mass (h-'M^) 



Figure 1. The virial mass function using cosmological simula- 
tions listed in Table [T] The error bars inside the symbols repre- 
sent the Poissonian error. The analytical mass function of Sheth 
and Tormen provides a very good fit for our data. 

and mass resolution (6.14 x 10^ H'^Mq, using 8.98M par- 
ticles) with the same cosmological parameters as the simu- 
lations mentioned above. It is clear that the density profile 
does not evolve in a si gnificant way from z = 1 to z = 0, as 
already pointed out bv lDiemand et al.1 (|2007l ). Nonetheless, 
the virial radius varies by about a factor of two as marked 
with vertical lines in Figure [5] This growth of the virial ra- 
dius is obviously an artifact of overdensity-based mass that 
does not reflect the real accretion of mass in the halo but 
the evolution of the mean matter density of the Universe, 

Pm(z) = (1 + z)^Pm,0- 

Therefore, a detailed analysis of the structure of dark 
matter haloes may prove to be valuable in the search for the 
virialized regions of the haloes. For example, the study of the 
radial phase-space diagram provides some insights on halo 
structure. In Figure |3] we show the radial velocity pattern of 
dark matter particles inside haloes of different masses taken 
from our simulation boxes. The top, middle, and bottom 
panels are for a low-mass halo with Mvir = 2.9x lO" /i'^Mq, 
a galactic-size halo of mass Mvir = 1.4 x 10^^/i~^Mq and a 
cluster-size halo with Afvir = 1.3x 10^^/i~^Mq, respectively. 
About 20 per cent randomly selected particles are shown 
for the low-mass halo and the galactic-size halo, and 100 
per cent for the cluster-size halo. Th e structure o f haloes in 
phase-space was already discussed in lBusha et al.l l|2005h . In 
the inner parts the average radial velocities are zero. Around 
the virial radius there are some signs of infall or outflow, 
as it will be discussed below. At large distances, the Hub- 
ble flow is the most remarkable feature, except for some 
spikes showing the presence of neighbouring haloes. This 
plot shows a substantial difference in the structure of dark 
matter haloes associated with different masses. Although the 
transition from the inner region to the Hubble flow shows 
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Table 1. The main parameters of the cosmological simulations. 



BoxSOS BoxSOG Boxl20 Box250c 



Box Size (h-'^Mpc) 




80 


80 


120 


250 


Mass of a particle (h~ 




4.91 X 10'' 


3.18 X 10** 


1.07 X 10^ 


9.67 X 10^ 


Spatial Resolution {h' 


-ikpc) 


0.52 


1.2 


1.8 


7.6 


Number of particles 




159.8M 


512^ 


512^ 


512^ 
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Figure 2. Evolution of the density profile for a galactic-size halo. 
The vertical lines mark the virial radius at each redshift. The 
change in the virial radius by a factor of two since redshift z = 1 
gives a false impression of a strong evolution of the halo. Yet, 
there was a very little evolution in the physical density in this 
halo since z = I. 

monotonically increasing radial velocity in the low-mass and 
the galactic-size halo, the situation is different for clusters 
where there is a significant tail in the velocity distribution 
with large negative velocities (see also Fig. 3 in lBusha et al.l 
This infall is altering the equilibrium of the particle 
distribution in the halo, and hence th at region is no longe r 
virialized, as already pointed out by iMaccio et al.l (|2003h . 
This suggests a natural choice for the boundary of a dark 
matter halo, which is defined by the innermost radius in 
which mean radial velocity is equal to zero. In practice, we 
use a small threshold in order to reduce effects of the statis- 
tical noise. The deviation of the spherically averaged radial 
velocity profile from zero is a signature of non-virialization. 
For this reason the zero- velocity radius is defining the mass, 
which is associated to a halo. This radius will be hereafter 
referred to as the static radius, and the mass inside a sphere 
of this radius, as the static mass of the halo. 

Although the threshold in the mean radial velocity is 
somewhat arbitrary and may affect our results, we take ad- 
vantage of the good statistics in our simulation boxes to 
reduce this threshold to only 5 per cent of the virial ve- 
locity Vvir = sj GMvii / Rvir of the halo. This allows us to 



get acceptable signal-to-noise ratios in the average profiles 
for all the mass ranges under study, without considerable 
overestimation of the virialized region. 

Individual haloes usually have similar radial profiles at 
a given range of mass. Thus, the analysis of the average halo 
profiles is reasonable. The radial profile of a given physical 
quantity is then obtained by getting the median value of this 
quantity over all the haloes in the sample for a given mass 
bin. This is done for every radial bin, which are logarith- 
mically separated for convenience. For the innermost radial 
bins, there is an additional constraint: the median profile 
is representative of the halo sample at a radial distance R 
only if a large fraction (about 80 per cent) of the haloes 
in th is sample contain at least 200 particles l|Klvpin et al.l 
bOOlT ) inside the sphere of radius R. 

As mentioned above, there is an important mass depen- 
dence in the radial velocity distribution of the haloes. The 
analysis of the median radial velocity profile reveals a promi- 
nent infall for cluster-size haloes and a noticeable outflow for 
low-mass haloes, whereas galaxy-size halo es show a balanced 
situation between infall and outflow (see IPrad a et aLll2006l 
for more details). This is represented in Figure |3] This gen- 
eral trend is dividing the total sample in two sets of haloes: 
those which are exceeding the threshold from below (infall) , 
and those which exceed it from above (outflow). 

Other properties of haloes also depend on halo mass. 
Figure [S] shows some of them for the same mass bins as in 
Figure m Upper panels represent the radia l and 3D veloc- 
ity dispersions (see also lWoitak et al]|2005l for an extensive 
analysis of radial v elocity moments in massive haloes, and 
IWoitak et al.ll200^ for a phenomenological model which de- 
scribes these). The radial velocity dispersion for haloes with 
higher concent ration has a larger inner value in terms of the 
virial velocity (|Lokas fc Mamonll200lf ). In the outer region 
there is a local minimum, which radius and magnitude also 
show dependence on halo mass. In the bottom left panel we 
show the median density profiles multiplied by (r/i?vir)^ in 
order to make more clear any change in their logarithmic 
slope. At smaller distances r < 0.1i?vir (not shown in the 
figure) curves start to decline because the slope is getting 
smaller than —2. It is also noticeable that beyond 1-2 i?vir 
the density profiles do not decrease as as for the NFW 
density profile. For a detailed parametrization far beyond 
7?vir see Tavfo et al. (in preparation). Finally, the bottom 
right panel shows the radial profile of the circular veloc- 
ity, which by definition coincides with the virial velocity 
at r = Rvir- The maximum of this profile is again mass- 
dependent, as it is somewhat related to the inner concentra- 
tion. 

The analysis of the mean radial velocity profiles over 
all masses (similar to those in Figure |4]) allows us to draw 
the relation between the static mass and the virial mass. 
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Figure 3. Phase-space diagram for tlie particles in dark matter 
haloes. The top panel represent a low-mass halo (Mvir = 2.9 X 
lO^^/i~^M0) and the central panel displays a galactic-size halo 
(Mvir = 1-4 X 10^'^/i~^Mq), whereas the bottom panel shows 
a cluster-size halo (Mvir = 1.3 x IO^^/i^^Mq). Clusters clearly 
show strong infall pattern around the virial radius, whereas this 
is not observed in smaller size haloes. 



0.8 



0.6 



> 0.4 



0.2 



-0.2 



10'=h-'Ma < M,,^ < 2*10'2h-iMs 
10"h-'M„ < M„,. < 2*10"h-'M„ 



/ ' 



/; 

// 

// 



+ 5% 



■■/■■;'■ 



-5% 



0.1 



10 



r/R , 



Figure 4. Mean radial velocity for three different mass bins. 
The profiles were obtained by averaging over hundreds of distinct 
haloes on each mass bin. In dotted line is shown the selected 
threshold delimiting the static region (5 per cent of the virial 
velocity). Cluster-size haloes display a region with strong infall 
(dashed line). On the contrary, low-mass haloes (solid line) and 
galactic haloes (long-dashed line) do not show infall at all but a 
small outflow preceding the Hubble flow. 



This relation is shown in Figure [6] Solid circles represent 
the Afstatic/Afvir ratio for the median profile built from all 
distinct haloes, while open circles show the same ratio for 
the median profile of isolated haloes. The isolation criteria 
here is that the nearest neighbour of a halo is further than 
three times the sum of the virial radii of both haloes. The 
uncertainty on this relation is difficult to estimate, as static 
masses could not be obtained for individual haloes. Our ra- 
dial velocity profiles are too noisy (because of infalling and 
outflowing subhaloes) to distinguish the static region, and 
hence we stacked these profiles in order to analyse them. 
Thus, we calculate the halo-to-halo variation of this rela- 
tion by measuring for each halo the mass inside the median 
static radius. In this way we obtain the Mstatic/Afvir ratio 
for isolated haloes, shown in this figure as dots. The distri- 
bution of this ratio has a long tail to large values of this ratio 
due to surrounding neighbours, which are inside the median 
static radius. Again, these unrealistic values of Mstatic/Afvir 
are due to the lack of a reliable determination of the static 
mass for individual haloes. There are no large differences 
between distinct and isolated haloes (except for those bins 
with poor statistics), so we take the Mstatic/Mvir relation for 
distinct haloes, as the isolation criteria reduces their number 
drastically. 

The dependence of this relation with the virial mass 
is straightforward. The outflow in low-mass haloes is 
less prominent with increasing mass, and hence the ratio 
Matatic/Mvir iucreasBS. A drop in this ratio occurs due to 
the transition to infall around 5 x 10^^/i~^Mq, and then it 
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Figure 5. Median profiles for the same halo mass bins as in 
Figure|4] These profiles show how the behaviour of haloes depends 
on the halo mass. Top left panel: radial velocity dispersion. Top 
right: 3D velocity dispersion. Bottom left: density profile. Bottom 
right: Circular velocity profile. The different line styles represent 
the same mass bins as in Figure l4l 



decreases due to higher infall with increasing virial mass. 
As it may be readily checked in Figure |4j the lowest mass 
in which infall takes place is very threshold-dependent. To 
alleviate this, we assume that this relation is a continuous 
function which is fairly well approximated by two power 
laws: 



logio(Mstatic/Mvir 



Qi logio(Afvir) + Pi region 1 
Q2 logj(,(Mvir) + /92 region 2 



(2) 



where the mass is measured in units of /i~^Mq. In Table [2] 
we list the values of the best-fitting parameters for a and /3. 

Interestingly, concentration of individual dark matter 
haloes decreases with static mass, in a similar way as it does 
with Afvir- This makes the ratio Mstatic/Afvir non sensitive 
to halo concentration. Again, the static mass of individual 
haloes is estimated by measuring the mass inside their me- 
dian static radius corresponding to its virial mass bin as in 
Figure [6l The Mstatic(c) dependence is similar to Mvir(c) 
because Mstatic is close to Mvir (the difference is less than a 
factor of two in the mass range under study). Hence corre- 
lations involving the logarithm of the mass can only change 
their zero point by ~ 0.3 if we use the static mass instead 
of the virial mass. The values of the parameter a in Table [2] 
suggest that the slope would present also a small change of 
~ 5% which might be important in some correlations. In the 
particular case of t he correlation betwee n halo mass and halo 
concentration (e.g. iMaccio et al]|2007l 'l. the replacement of 
virial mass to static mass introduces a modification which is 
small enough so that this correlation remains almost entirely 
unaffected. 
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Figure 6. The halo static mass - virial mass relation at z = 0. 
Dots show the relation for individual isolated haloes. Circles show 
this relation for isolated (open) and distinct (solid) median haloes 
at different mass bins. The ratio Matatic/A^vi^ increases due to 
reduced outflow up to the first appearance of infall, making this 
ratio to decrease. Dashed line shows the approximation to two 
power laws connected where they cross each other, the parameters 
are those listed in Table |2] 



4 ANALYSIS OF EQUILIBRIUM 

Haloes of different mass clearly show very different structure. 
At small (r < R^ir) distances the main differences are due to 
the differences in concentration. The haloes inside the virial 
radius are nearly in equilibrium (with possible exception of 
the most massive clusters). At larger distances the differ- 
ences are more com plicated. Figure |4] and similar results in 
iPrada et al.l (|2006l ) demonstrate that haloes of small mass 
have static radius few times larger than the formal virial 
radius. To make things even more complicated, the small 
haloes do not have an infall region. We can try to somewhat 
clarify the situation by applying the Jeans equation. In order 
to do this, we first need to measure the velocity anisotropy 

In Figure [7] we show the median radial profile for {3 for 
the same mass bins as in Figure |4] Inside the virial radius 13 
is positive. Cluster-mass haloes have a familiar behaviour: (3 
increases with distance and gets to /? ~ 0.3 — 0.4 around the 
virial radius. The real surprise came for the small galaxy-size 
haloes. After initial increase, (3 reaches maximum (3 « 0.2 
and then declines with distance. It goes almost to zero at 
r = (1 — 2)i?vir. At larger distances the velocity anisotropy 
is negative and approaches zero at very large distances. We 
should note that (3 = does not mean nearly circular orbits. 
Orbits of dark matter particles or satellites are typically very 
elongated even for isotropic velocity ellipsoid: axial ratios are 
typically 1:4-1:5 for /? = and the NFW profile. 

Why small haloes do not have a turn-around radius fol- 
lowed by an infall region? The top-hat model of non-linear 
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Table 2. Parameters of the approximation of the Mgtatic/A^vir ratio as a function of the virial mass. 



region 




mass range 




a 




/3 


1 
2 


1.0 X 10^0 
7.4 X IQii 


< Afvir//i-^M0 


^ 7.4 X 10" 
^ 5.0 X 10" 


oi = +0.052 ± 0.005 
Q2 = -0.055 ± 0.009 


/3i = 
& = 


-0.40 ±0.05 
+0.87 ±0.12 



0.5 

0.4 
0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 
-0.4 



-0.5 




10'=h-'M„ < M,.,. < 2»10i2h-'M, 



10"h-'M3 < M,,^ < 2*10"h-'Mo \ / 



0.1 



10 



r/R„ 



Figure 7. Velocity anisotropy profile for the same mass bins as in 
Figure |4] and Figure [5] In the inner region the velocity anisotropy 
is preferentially radial. It changes to nearly isotropic in the outer 
regions with (1 — 2)iJvir followed by preferentially tangential 
at larger distances. 



evolution is supposed to provide a reasonable description of 
what happens at large distances. Yet, it seems to be failing 
qualitatively for galaxy-size haloes, which do not have infall 
region and even show an outflow. Haloes are complex sys- 
tems: they are not spherically symmetric and may have sub- 
stantial substructure. Still, the spherically symmetric non- 
stationary Jeans equation is an undeniably significant im- 
provement as compared with the top-hat model. The Jeans 
equation provides a simple approximation for the haloes. As- 
suming spherical symmetry the non-stationary Jeans equa- 
tion can be written in the following form: 



r d{pVr 
P 



dt 



+ (Tr 



dlnr 



+ 2/3 



= Vc (r). 



(3) 



where Vr is the radial velocity and v^{r) = GM{< r)/r. 
When the rms velocities are negligible, the equation Q 
is the normal equation of motion. Assuming a small initial 
perturbation and the Hubble flow, we recover the top-hat 
model. Yet, the haloes do not have small rms velocities even 
well outside the formal virial radius. Figure |8] illustrates the 
point. Even for the massive haloes, which typically dom- 
inate their environment, the rms velocities are relatively 
large. The effect is even stronger for less massive haloes. 
Another important component is the density gradient. Fig- 



ure [2] shows that the density gradients exist at very large 
distances: at many virial radii. At z = 2 the density was on 
average steeply declining well outside its virial radius. 

The left-hand-side of the Jeans equation has two terms. 
The first term describes non-stationary effects. It is the dom- 
inant term at early stages of evolution. The second term is 
related with the pressure gradient. It dominates in the virial 
part of the halo. If the system is stationary, the second term 
must be equal to the square of the circular velocity in the 
right-hand-side of the Jeans equation. We can call this term 
the equilibrium circular velocity: 



null 



dln{pa^ 



din? 



+ 2/3 



(4) 



We plot Vc.equii in Figure [8] For galaxy-size haloes (right 
panel) in the inner region r < (2 — 2.5)i?vir the term is very 
close to the right-hand-side of the Jeans equation, which 
means that the system is nearly in virial equilibrium. At 
larger distances the term falls short of the circular velocity: 
the halo is out of equilibrium. 

For clusters (the left panel) the equilibrium is only in 
the central r < R^ir region. Just outside of this region at 
r = (1 — 2)i?vir the equilibrium term Uccquii is substantially 
above the circular velocity Vc- At the same range of radii the 
infall velocity is clearly detected. The upturn in the Hccquii 
is the cumulative result of two effects: increase in the slope 
of the density profile and substantial change in the velocity 
anisotropy. The same trend exists in the case of small haloes, 
but it is barely visible. At larger distances r > 2iivir the 
equilibrium term Uc.cquii falls below the circular velocity. 

Our analysis of the Jeans equation shows that the mean 
radial velocity is a good indicator of the virialization of a 
halo and that the low-mass haloes are nearly in virial equi- 
librium well beyond formal J?vir. 



5 EVOLUTION AND SCALING OF HALO 
MASSES 

Dark matter haloes grow by accreting mass from their out- 
skirts. Haloes with small mass grow very slow at late times. 
Thus, in the past they must have experienced a remark- 
able infall, which has created extended virialized regions sur- 
rounding the haloes at present. In this section we study how 
this process of building of halo exterior was happening. For 
this reason we analyze our simulations at different redshifts. 
We then take the median of different halo properties using 
the same mass bins as in previous sections with masses rang- 
ing from 10^° to 10^^/i~^Mvir. As expected, the statistics in 
each mass bin is smaller at higher redshifts. Using these data 
we calculate the relation between static and virial mass for 
z = — 2. When doing this analysis we take every halo at 
a given redshift without considering their relation with the 
final haloes at z = 0. The result is shown in Figure |9l where 
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Figure 8. Different velocity components for cluster-size haloes 
(left panel Mvir si 2 X 10^"^ h~^yiQ) and galaxy-size haloes (right 
panel M^ir ~ 10^^/i~^Mq). The dot-dashed curves shows the 
predictions fcequil of stationary Jeans equation Q. The sta- 
tionary solution closely follows the real circular velocity up to 
(2 — 2.5)-Rvir for galaxy-size haloes. It falls below Vc at larger 
distances indicating significant non-stationary effects. The situa- 
tion with clusters is different: the central virialized region is sur- 
rounded by a shell where virialization is happening {vc^equH > f c) 
followed by the region where "Uceijiiii < ""c- See text for details. 

we plot only the mass bins with enough statistics. At z = 
the Mstatic/Mvir ratio increases with mass. Then it reaches 
mciximum and starts to decline at large masses. Because of 
mass resolution we do not see the rising part of Afstatic/Mvir 
ratio at large redshifts, but the declining branch of the curve 
is better resolved, and we clearly see that the ratio gets be- 
low unity at large masses. 

It is suggesting that the relation between static and 
virial mass is essentially scale-free for different redshifts up 
to z = 2. In order to test this hypothesis, we approximate the 
declining part of this relation by a power law Afstatic/Mvir = 
CM"i^, and then find the crossing point: Mstatic ~ M^ir- 
This solution Aio = Mo{z) allows us to scale the static to 
virial mass relation for different redshifts. The results clearly 
overlap, indicating that there a universal (z— independent) 
static to virial mass relation shown in Figure [TOl This plot 
shows that the static to virial mass ratio is well fitted by 
two power-law with different slopes: 

M,tatic _ / lO+O-^^^iO O^l ^-H0.054±0.005 ^ < ^g-S 
1Q-0.010±0.005 ^-0.075±0.004 > IQ-^ 

where x = Mvir /Mo . This relation allows us to convert the 
virial mass, which is easier to be computed accurately in 
numerical simulations, to the static mass. In order to have 
a complete description of this relation for any redshift, it 
remains to be determined the redshift dependence of the 
mass scale Mo{z). Although this dependence may suffer from 



Figure 9. The static to virial mass relation from 2 = 2 to 2 = 0, 
taking into account every halo at a given redshift. Large haloes 
having Mgtatic < ^vii exist at all redshifts, but the transition 
mass decreases with increasing redshift. 




10-e 10-5 10-4 10-3 10-2 10-1 10° 10' 10= 10^ 



Figure 10. Same as Figure|9l but scaled to the mass Mo in which 
the fit of the smooth decline in Mgtatic/Afvir to ^ power law is 
equal to one. The time dependence of the Mstatic^Afvir relation 
is encoded in the mass scale A/q. The solid line represents a fit 
to a power law at both sides of the transition mass, see equation 
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Figure 11. The evolution of the mass scale Mq. This behaviour is 
approximately fitted by a power-law, Mq ~ 2.8 X lO^^a*'^/i~^M0 
(see text for details). 



the scatter in the determination of the static mass, it seems 
clear from Figure[TT]that there is an approximate power-law 
dependence Afo(a) = Mo{a = l)a' , from z = 2 to z = 0. 
The best fit to our dat4i] is given by: 



Mo 



(2.8 ±0.8) X lO^^a** 



*h-'M 







(6) 



It is worth of pointing that the absence of explicit red- 
shift dependence in the relation between static and virial 
mass when the variable M^ii/Mo is used, allows an ana- 
lytic and straightforward derivation of the relation between 
both mass accretions, i.e. taking the logarithmic derivative 
in equation (O and using equation ^ we obtain: 



Mstatic ^ " Mvir a ' 



(7) 



where the dot refers to the derivative with respect to the 
scale factor a. However, this equation is no longer valid near 
the transition mass, where this derivative is not finite. 

Despite the mass scale Mq is not motivated by the 
physics of the formation of dark matter haloes (as opposed 
to e.g. M,), it is the best way we found to encode the redshift 
dependence of A/static /Mvir. We have also explored the rela- 
tion of this ratio with v = 5c/cr{M, z). However, we found a 
wider scatter for a given value of v, with lower values of the 
ratio for increasing redshift, in particular for v > 1. 



^ We want to emphasize that halo virial mass is defined here 
with a overdensity criterion which evolves with time, as detailed 
in Section [2] The choice of non-evolving overdensity criterion will 
imply a change in the Aio{z) relation. 



6 MASS FUNCTION 

The relation between the static and the virial mass is not 
monotonic, as we have shown before. This implies a change 
in the shape and in the amplitude of the distribution of 
the number density of haloes with a given mass, i.e. the 
mass function. As the virial mass is calculated accurately 
for all the haloes in our simulations, we may also obtain 
their static mass directly from the conversion formula given 
in equation ([Sjl. This allows us to estimate the static mass 
function of dark matter haloes and compare it to the virial 
mass function. We obtain very similar results taking the 
virial mass bins in Figure [10] and then shifting them accord- 
ing to their corresponding value of the ratio Mstatic/Mvir- 
The results are presented in Figure 1121 In spite of the fact 
that the mass function of ISheth fc TormenI (1999) provides 
a very good fit for the virial mass (as shown in Figure [T] for 
the case of z = 0), we find that this analytic model is not 
a good description for the static mass function. Moreover, 
some authors using different definitions of virialized halo 
have already shown the presence of significa nt departures 
from the Sheth & Tormen function, see e.g. Im accio et al.l 
(l2003l). 



Surprisingly, the mass function of iPress fc Schecht"eil 

1 19741 ) approximates really well our data at « = in the 
mass range under study. In order to test robustness of our 
estimates of the static mass function, we changed different 
parameters (e.g., changed the threshold to O.lVvir in algo- 
rithm of detecting the static radius). Results are stable. We 
leave for future the inves tigation of the static mass func- 
tion at its high mass end l|Betancort-Riio fc Montero-Dorta 
|2006| ). because the poor statistics for the highest mass bins 
here prevents us from drawing any conclusion about it. The 
redshift evolution of the mass function is also shown in Fig- 
ure 1121 We realize that the static mass function deviates 
from the corresponding Press & Schechter function ai z = 1 
and z = 2 so it is no longer described by this model. None 
the less, it is interesting that we can still recover the agree- 
ment with Press & Schechter by an overall increase in the 
static mass. This shift could be accounted for by the scatter 
in the Mq (z) relation, which is very similar to this deviation. 
However it is not clear whether this deviation is an artifact 
due to resolution effects, so this remains as an open issue. 



7 EVOLUTION OF THE MAJOR HALO 
PROGENITOR 

We are now in a position to study the evolution of the mass, 
i.e. the mass accretion history, of dark matter haloes in cos- 
mological simulations. So, we need to track the set of haloes 
identified by EDM at z = back in time and then measure 
their static mas s . Her e we use a method similar to that in 
lAscasibar et al.l l|2007l ). First of all, we pick up a sample of 
few hundreds of particles around the density maximum of 
every halo at z; = 0. These particles are then tracked back 
to the highest redshift available in our cosmological simula- 
tion. Although it may occur that these particles no longer 
belong to the same object at that redshift, we take the par- 
ticle which has a minimum sum of squared distances to the 
others in the sample, as an initial guess for the centre of the 
protohalo. Finally, we calculate the centre of mass around 
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Figure 12. The static mass function obtained from individual 
haloes using equation ^ (solid circles for 2 = 0, triangles for 
2 = 1 and squares for z = 2). Error bars show the Poissonian 
error only (not including the uncertainty in the Mgtatic ~ ^vir 
relation). Solid and dashed lines show the PS and ST mass func- 
tion, respectively. Surprisingly, the PS mass function provides a 
very good fit to our data at 2 = 0. On the other hand, the mass 
function at higher redshifts deviates from the PS analytic model. 



this position iteratively using the whole set of particles in 
the simulation box. After few iterations, we obtain the cen- 
tre of the primordial peak of density which corresponds to 
our halo at z = 0. We then label some particles around this 
centre. This new set of particles will help us, by calculation 
of their centre of mass, to identify at every redshift the po- 
sition of the halo which encloses them. Hence this halo, as 
the offspring of the main primordial peak corresponding to 
the halo at present, will be referred to as the halo major 
progenitor. 

In our analysis we study the evolution of the pro- 
genitors from three different mass bins at z = 0: 
low-mass haloes (lO^'^"-lO^'^'^/i~^M0), galactic-size haloes 
(lO^^°-lO^^-^/i"^M0) and cluster-size haloes (lO^^'^- 
10^'*"/i~^Mq). These ranges are selected in order to study 
the dependence of mass on our results with very good statis- 
tics: each bin contains more than a thousand of haloes. This 
allows us to calculate with good accuracy the median pro- 
files for different physical quantities at every redshift from 
z = to 2 ~ 2. We keep this as the maximum redshift anal- 
ysed here because the median radial profile is no longer well 
resolved (i.e. with more than 200 particles for the major- 
ity of the haloes) in the inner regions for higher redshifts. 
Thus, we cannot determine the size of the virialized region 
and hence we cannot go further in redshift. 

The evolution of the mean radial velocity profile for 
the median halo is displayed in Fi gure fT3l Here we s how 
that the time sequence suggested in iBusha et all (|2005l ) for 
cluster-size haloes takes place even for low-mass haloes: the 



Figure 14. The static to virial radius as a function of the redshift. 
It appears that this ratio behaves differently for each mass bin. 
The mass bins are the same as those on Figure 1131 Low-mass 
haloes have reduced their static region in units of i?vir due to 
recent outflow. On the other hand, galaxy-size haloes are reaching 
a maximum on this ratio. Cluster-size haloes continue increasing 
this ratio as the infall is less prominent. 

mean radial velocity profile experiences infall in their outer 
parts while forming. Just after the infall is finished, there 
is an epoch in which some outflow appears, followed by a 
final profile with an inner virialized region and the Hubble 
Moreover, the time scale o f this sequence see ms to 
be mass-dependent: the analysis of iBusha et al.l (|2005l ) finds 
this sequence in the future for cluster-size haloes, whereas 
low-mass haloes are already in the outflow phase at present. 
Galaxy-size haloes are now experiencing neither infall nor 
outflow, so the virialized region is reaching a maximum in 
units of the virial radius, as already shown in Figure [S] 

This evolution is reflected in the redshift dependence 
of the static radius. In Figure [14] we illustrate the evolution 
of the ratio of static to virial radius. It is clear from this 
plot that low-mass haloes have already formed (the infall 
stopped when the maximum ratio took place) and they are 
even losing some of their virialized mass at present. Galaxy- 
size haloes are reaching now this maximum, but for cluster- 
size haloes the infall has not even flnished. 

The static mass accretion history is presented in Fig- 
ure [TS] The solid, long-dashed and short-dashed lines show 
the evolution of the static mass of the median halo for three 
different mass bins, while the dotted line corresponds to 
their virial mass. The error bars are estimated using the 

^ This final state agrees with the picture shown by other dif- 
ferent analyses of halo virialization, in which virial parameters 
are ill-defined in merging events but there is a general trend for 
the virial ratio 2Ei;in/|Epot | to approach unity at present, e.g. 
iHetznecker fc Burkertl booel) . 




Figure 13. The evolution of the mean radial velocity profile from z = 2 to z = 0. Left panel: low-mass haloes. Middle panel: galactic-size 
haloes. Right panel: cluster-size haloes. The sequence followed by low-mass haloes has hardly started by cluster-size haloes. 



halo-to-halo variation, calculating the mass inside the me- 
dian virial radius for each dark matter halo, and taking the 
central 68 per cent of these values. Neither the uncertainty 
in the virial radius nor the scatter in the static vs. virial 
mass relation were taken into account. Despite of the scatter, 
some conclusions may be drawn from this plot. The recent 
assembly history is very mass-dependent: low mass haloes 
hardly experienced a change in their mass since z = 1, and 
their virial mass is under-predicting the average amount of 
matter inside them since z ~ 2. On the other hand, galactic- 
size haloes have been increasing their static mass up to the 
present epoch, most likely due to halo relaxation. This pro- 
cess is responsible of the largest difference between the static 
and the virial mass at 2 = 0: the mass inside the virialized 
region is about a factor of two larger than the virial mass. 
Finally, cluster-size haloes are still increasing their mass at 
present, but the difference with respect to the virial mass is 
not so large. In fact, around z — 0.5 both static and virial 
mass should have been very similar. On the contrary, it is 
likely that they differ significantly at higher redshift. 

We find for our me dian profiles th a t the function pro- 
posed in equation (3) in lWechsler et akl (|2002h for the virial 
mass, Mvir(a) ~ Mvir(a = 1) exp {—5z) provides an excellent 
fit for the accretion history of the virial mass from 2 = 2 to 
2 = 0. However, if we include this relation in our model 
given by equation (0, and solving the differential equation 
we obtain: 



MBtatic(a) 
Mstatic(a = 1) 



exp 



-(l + a)Sz 



(8) 



which fails to reproduce the evolution of the static mass 
when the discontinuity of the static vs virial mass relation 
takes place. This is reasonable as equation ((Tj is obtained 
when all the haloes (not just the progenitors of the haloes 
at 2 = 0) are taken into account. None the less, this starting 
point suggests the functional dependence of this evolution, 
so that a fit of the data to: 



Mstatic(2) = Mstatic(2 = 0)(1 + 2)" 



(9) 



where both a and /? are free parameters, provides a better 
description of the data. This function is very similar to that 
corresponding to the virial mass, but there is an extra factor 
which enhances the fit, in particular for low-mass haloes. 
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Figure 15. The mass accretion history of dark matter haloes. 
The evolution of the static mass is indicated by the crosses 
and the error bars, with their corresponding fit to equation (|9]l 
(solid, long-dashed and dashed lines for low-mass, galaxy-size and 
cluster-size haloes, respectively), while the dotted line shows the 
evolution of the virial mass, according to the parametrization by 
IWechsler et aL I ||2002D . The mass bins are the same as those on 
Figure [13] The static mass inside low-mass haloes is nearly con- 
stant from z = 0.5 to z = 0. On the other hand, galaxy-size haloes 
increase their static mass even steeper than their virial mass, due 
to relaxation around haloes. Clusters are accreting mass, but also 
incorporating some of the recently virialized mass from their sur- 
roundings. 



As we pointed out above, the recent increase in the mass 
of galactic-size haloes is mainly due to virialization of mass 
surrounding them. In Figure [16] we display the evolution of 
the density profile for the median galactic-size halo, which 
is similar to the profile for individual haloes as shown in 
Figure [2] Only small variations in the density profile occur 
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Figure 16. The evolution of the median density profile from 
z = 2 to z = for galactic-size haloes. As shown in Fig. |2] for 
individual haloes, there has been very little change in this profile 
from z = 1 up to the present epoch. The vertical dotted lines 
marks the virial radii at 2 = 2 (left), z = I (middle) and z = 
(right). Static radii are shown with vertical solid lines. The growth 
of the static radius in Milky- Way size haloes from 2: = 1 to 2 = 
is not due to accretion but virialization of the surrounding regions 
of the halo. 



for the median galactic halo from z — 1 to z — 0, which 
correspond to a small variation of the mass inside a phys- 
ical radius. However, these haloes experience an important 
growth of their static mass due to the expansion of the viri- 
alized region, as it is clear from Figure [13] 



8 DISCUSSION 

A more realistic (physical) estimate of the halo mass might 
be crucial for some areas in cosmology which are sensitive 
to accurate values of the mass of dark matter haloes. For 
example, our results may be relevant for understanding the 
galaxy formation processes that happen inside dark matter 
halo es. Indeed, semi-ana lytical models of galaxy formation 
(e.g. ICroton et aP l2006h might need to assess the impact 
of their predictions when the virial mass is no longer used. 
If the static mass and its evolution are used, one may get 
a more realistic picture of the structure formation process 
drawn from cosmological simulations. Moreover, recipes for 
star formation, AGN and supernova feedback and merging 
of galaxies are often taken from quantities of their corre- 
sponding dark matter haloes, such as their mass. Thus, the 
conclusions (e.g. colour-magnitude diagram, evolution of the 
stellar mass, etc.) from the analysis of these simulations us- 
ing semi-analytical models of galaxy formation may be bi- 
ased due to an underestimation of the mass associated to 
a halo if the former virial mass is considered, especially for 



the galaxy-size haloes such as the host of the Milky Way as 
we have reported in this work. 

The precise determination of the halo mass function, in 
terms of the static mass, can also be very relevant. Mass 
function of cluster-size haloes must be predicted with spe- 
cial accuracy because the number density of galaxy clus- 
ters is a probe of fundamental cosmological parameters, such 
as dark energy. This is very important since the evolution 
of the mass function of galaxy clusters is taken as an al- 
most direct measu rement of the equation of state of dark 
energy ui = P/p (|The Dark Energy Survey Collaboration! 
I2OO5I) and any statistically significant deviations in the mea- 
surements of the mass function of clusters from the predic- 
tions of the ACDM model, may be and likely will be in- 
terpreted as deviations of the dark energy from a simple 
cosmological constant. 

Other possible observational implications of this work 
also deserve further investigation. It is noteworthy that the 
transition mass of ~ 5 x 10'^^/i~^Mq, that separates the in- 
fall and outflow reg i mes, i s similar to the mass scale fo und by 



tall and outHow reg i mes, i s similar to the mass scale lo una by 
iDekel fc BirnboimI (|2006l ). iDekel fc BirnboimI (|2006l ) argue 
that below M < IO^'^Mq the efficient early star formation 
may occur due to cold fiows, as thermal pressure is not large 
enough to avoid the gas accretion from the outer gas reser- 
voir to the inner core of the halo. On larger scales the heating 
due to virial shocks may shutdown the star formation. This 
different behaviour of the gas around dark matter haloes in 
si mulations above and below this mass scale was proposed 
bv IDekel fc BirnboimI ^M) 

as an explanation for the bi- 
modality observed in many properties of galaxies such as 
the bimodality seen in the colour distribution of SDSS field 
galaxies. Thus, a given velocity pattern of infall or outfiow 
in a dark matter halo of mass above and below 10^^ h~^MQ 
may be correlated with the fact that the galaxy within the 
halo belongs to the blue cloud or the red colour sequence, 
respectively. Yet, this is only an speculation which remains 
to be investigated. However, the mass scale for transition 
between infall and outfiow is the same as the mass scale at 
which baryonic matter cannot penetrate the inner parts of 
the halo. At first glance this is just a matter of coincidence, 
as for higher redshift s the transition mass wou ld be well be- 
low the mass scale o f iDekel fc BirnboimI l|2006l ). Still, it is an 
interesting question whether the haloes experiencing infall 
of the dark matter also host galaxies, which grow fast by 
accreting baryons in their central regions. 

Such correlations between galaxies and their host haloes 
could unveil some of the relations between dark and baryonic 
matter, and hence cast some hght in the evolution of galaxies 
in the expanding Universe. 



9 CONCLUSIONS 

In this paper we use a set of high resolution A'^-body cos- 
mological simulations for the analysis of the mass inside the 
region of dark matter haloes with no net infall or outfiow 
velocities, i.e. the static region. Haloes show a typical ra- 
dial velocity pattern which depends on their mass: low-mass 
haloes tend to display a region with outfiow in their sur- 
roundings, while cluster-size haloes show a prominent infall 
velocity pattern. Galaxy-size haloes show an intermediate 
pattern with a sharper transition between the static region 
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and the Hubble flow. We find that the former virial radius 
tends to underestimate the size of the region with zero mean 
radial velocity for haloes with masses 

10^*/i~^Mq at 2 = 0. The mass inside this static region can 
be about a factor of two larger than the former virial mass, 
when the threshold defining the static radius is 5 per cent 
of the virial velocity Vvir = y^GMvir/ Rvii)- Lower values of 
the threshold (< 1 per cent) may lower this factor. However, 
we observe a clear trend in the Mstatic/Mvir vs A/vir relation 
which remains unchanged with this threshold: for low-mass 
haloes, this ratio is an increasing function of Mvir due to re- 
duction in the outfiow. On the other hand, high-mass haloes 
show a larger infall for larger Mvir, making this ratio to de- 
crease. The maximum occurs at Mvir — 5 x 10^^^~^Mq, 
where the size of the static region is the largest in units of 

-Rvir. 

The mass function of objects whose mass is defined in 
the way presented here resembles that of Press & Schechter 
in the range of mass we have studied, but only at z = 0. 
At higher redshifts the static mass function deviates signif- 
icantly from it. This disagreement could be accounted for 
by the scatter in Mo{z), although it seems more likely that 
the resemblance with Press & Schechter at 2; = is due to 
coincidence (in both shape and amplitude to within ~ 20%) 
of the ratio fps/fsT with Matatic/Mvir. Whatever the case 
may be, this resemblance is theoretically unexpected be- 
cause the Sheth & Tormen function, which fits very well 
the meiss function for virial masses, is derived from a more 
realistic approach (ellipsoidal collapse) . In any case it seems 
clear that when the virial mass is used, the number density 
of dark matter haloes of a given mass at z = might be 
underestimated in the mass range 10^''^"^Mq < Mvir < 

The redshift evolution of the Mstatic-Mvir relation turns 
out to be very weakly redshift-dependent from 2 = 2 to 2 = 
when appropriate variables are used. This dependence on the 
redshift is encoded as an evolving mass scale Mo which indi- 
cates the approximate mass at which Afstatic = Afvir in the 
declining part of the relation. By doing this, it is straight- 
forward to derive from this relation a differential equation 
for the evolution of the static mass, provided that the evo- 
lution of Mvir and of Mo are already known. As we have 
shown in this paper, the evolution of the mass scale has an 
acceptable fit to a power law: Mo ~ 2.8 x W^'^a^-^h'^MQ 
. On the other hand, the evolution of the virial mass Mvir 
of the major progenitor is w ell described by an ex ponen- 
tial law in 2, as claimed by IWechsler et al] (|2002l ). Once 
we have fixed the evolution of Mvir and of Mo, we solve 
this differential equation and we find a simple model for 
the evolution of the static mass of the major progenitor, i.e. 
Mstatic(2) = A/static(2 = 0) ( 1 -|- 2) "''e""^ , which turns out to 
be a very good fit to our data. In any case, there is a differ- 
ent behaviour for the evolution of the static mass in differ- 
ent mass bins. The static mass of low mass haloes is nearly 
constant from 2 — 0.5 to 2 = 0. Galaxy-size haloes keep 
growing at present, but the reason here is not the decreas- 
ing background density of the expanding Universe. Instead, 
the relaxation of the matter in the surroundings of these 
haloes is incorporating mass at even a higher rate than the 
unphysical growth of the virial mass. The same stands for 
clusters, although the static and virial mass are not very 
different at present. 



The mean radial velocity profile (averaged over a h alo 
mass bin) evolves in the way found bv lBusha etHI (|2005l ) for 
cluster-size haloes: while the halo is forming, it accretes mass 
from its surroundings via infall, which is being reduced with 
time. This is followed by a period of outflow of unbound mat- 
ter, and a final configuration thereafter consisting in a static 
region whose size is independent of time. Although the sim- 
ulation used in that paper, was evolved until a — 100, and 
our halo sample was instead tracked up to redshift 2 = 0, we 
find a situation that is consistent with this picture. More- 
over, we find that the timescale of this sequence is different 
for each mass bin: low-mass haloes at present are in the 
outfiow phase, but cluster-size haloes are still in their infall 
phase. Interestingly, galactic-size haloes are experiencing a 
transition (infall to outflow) epoch at 2 = 0. 
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